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Control over all internal and external degrees of freedom of molecules at the level of sin- 
gle quantum states will enable a series of fundamental studies in physics and chemistry 1 ' 2 . 
In particular, samples of ground-state molecules at ultralow temperatures and high number 
densities will allow novel quantum-gas studies 3 and future applications in quantum infor- 
mation science 4 . However, high phase-space densities for molecular samples are not readily 
attainable as efficient cooling techniques such as laser cooling are lacking. Here we pro- 
duce an ultracold and dense sample of molecules in a single hyperfine level of the rovibronic 
ground state with each molecule individually trapped in the motional ground state of an 
optical lattice well. Starting from a zero-temperature atomic Mott-insulator state 5 with opti- 
mized double-site occupancy 6 , weakly-bound dimer molecules are efficiently associated on a 
Feshbach resonance 7 and subsequently transferred to the rovibronic ground state by a stimu- 
lated four- photon process with >50% efficiency. The molecules are trapped in the lattice and 
have a lifetime of 8 s. Our results present a crucial step towards Bose-Einstein condensation 
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of ground-state molecules and, when suitably generalized to polar heteronuclear molecules, 
the realization of dipolar quantum-gas phases in optical lattices 8-10 . 

Recent years have seen spectacular advances in the field of atomic quantum gases. Ul- 
tracold atomic samples have been loaded into optical lattice potentials, allowing the realization 
of strongly-correlated many-body systems and enabling the direct observation of quantum phase 
transitions with full control over the entire parameter space 5 . Molecules with their increased com- 
plexity are expected to play a crucial role in future generation quantum gas studies. For example, 
the long-range dipole-dipole force between polar molecules gives rise to nearest-neighbour and 
next-nearest-neighbour interaction terms in the extended Bose-Hubbard Hamiltonian and should 
thus lead to novel many-body states in optical lattices in the form of striped, checkerboard, and 
supersolid phases 8-10 . 

An important prerequisite for all proposed molecular quantum gas experiments is the capa- 
bility to fully control all internal and external quantum degrees of freedom of the molecules. For 
radiative and collisional stability, the molecules need to be prepared in their rovibronic ground 
state, i.e. the lowest vibrational and rotational level of the lowest electronic state, and prefer- 
ably in its energetically lowest hyperfine sublevel. As a starting point for the realization of novel 
quantum phases, the molecular ensemble should be in the ground state of the many-body system. 
Such state control is only possible at ultralow temperatures and sufficiently high particle densi- 
ties. While versatile non-optical cooling and slowing techniques have recently been developed 
for molecular ensembles 11 , the achievable molecular phase-space densities are still far away from 
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the point of quantum degeneracy. Here, we exploit the fact that high phase- space densities can 
readily be achieved for atoms and that atoms can efficiently be associated on Feshbach resonances 
to form molecules 7 with minimal loss of phase-space density when an optical lattice is present. 
Subsequent state transfer to a specific hyperfine sublevel of the rovibronic ground state by means 
of a stimulated multi-photon process then preserves phase-space density and hence the quantum- 
gas character of the molecular ensemble. This approach is expected to allow the preparation of a 
molecular ground- state BEC 12 . 

A crucial ingredient for our experiments is the presence of an optical lattice. It provides 
full control over the motional wave function and prevents collisional loss. It allows us in par- 
ticular to maximize the efficiency for initial molecule production and the efficiency for ground- 
state transfer. In addition, our preparation procedure directly leads to a lattice-trapped molecular 
many-body state, which forms the basis for the many-body and quantum information experiments 
envisioned 810 . In the quantum gas regime without the use of an optical lattice, molecular state 
transfer to deeply-bound rovibrational levels of the singlet X S ground-state potential has recently 
been implemented for Cs 2 13 and KRb 14 . For KRb, the rovibronic ground state was reached, result- 
ing in a near-quantum-degenerate gas of fermionic ground-state molecules 14 . 

Our molecular quantum-gas preparation procedure is summarized in Fig. [T] We load a BEC 
of Cs atoms 15 into a three-dimensional optical lattice and drive the superfluid-to-Mott-insulator 
phase transition 5 under conditions that yield the highest number of lattice sites at which there are 
precisely two atoms (see the Methods section). We thus aim to maximize the size of the two-atom 
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Mott shell 5 . With up to 45% of the atoms at doubly-occupied lattice sites we come close to the 
theoretical limit of 53% given the parabolic density profile of the BEC 6 . The atom pairs reside 
in the motional ground state at each well and are then associated with near-unity probability 16 to 
Cs 2 Feshbach molecules, which are subsequently transferred to the weakly-bound level |1 >, the 
starting level for the optical transfer (see the Methods section) 131718 . Atoms at singly-occupied 
sites are removed by a combination of microwave and optical excitation 16 . We now have a pure 
molecular sample with near-unity occupation in the central region of the lattice. Each molecule is 
in the motional ground state of its respective well and perfectly shielded from collisional loss. 

We employ stimulated Raman adiabatic passage (STIRAP) 19 involving four laser transitions 
to coherently transfer the molecules into the lowest rovibrational level |5>= |t> = 0,J = 0>of 
the ground state singlet AT X £+ potential as shown in Fig. [2^, bridging a binding energy of hex 
3628.7 cm -1 « h x 109 THz 13 . Here, v and J are the vibrational and rotational quantum numbers, 
respectively, h is Planck's constant and c is the speed of light. For Cs 2 , a homonuclear molecule, the 
four-photon process is preferred to a direct two-photon process because it allows us to overcome 
small Franck-Condon overlap. Lasers L\ through L 4 couple |1 > and |5 > via three intermediate 
levels |2 >,|3 >, and |4 > (see the Methods section). For STIRAP in the presence of the lattice, 
the lattice light must not impede the transfer through optical excitation or by creating unwanted 
coherences. Also, the lattice wavelength has to be chosen such that the dynamical polarizabilities 
for 1 1 > and 1 5 > are closely matched in order to avoid excitation into higher motional states of 
the lattice as a result of motional wave-function mismatch 20 . We typically set the lattice depth 
to a value of 20 E R for atoms, corresponding to 80 E R for Feshbach molecules with twice the 
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polarizability and double the mass and 83 E R for molecules in v = at a lattice wavelength of 
1064.5 nm, as determined below. Here, Er (Er) is the atomic (molecular) recoil energy. 

Our experimental configuration ensures that only one particular molecular hyperfine sublevel 
is populated. The atomic BEC is prepared in the lowest hyperfine sublevel \F a = 3, m Fa = 3 >, 
where F a and m Fa are the total atomic angular momentum and its projection on the magnetic field. 
Feshbach association and transfer between Feshbach levels via avoided crossings, as illustrated 
in Fig. (see the Methods section), conserve 7 the total angular momentum projection M F = 
ui Fai + m F a 2 — 6. Fig. shows the hyperfine structure of the target state, i.e. the rovibronic 
ground state X l Y^ \v = 0, J = >. It splits into 28 hyperfine sublevels in the presence of a weak 
magnetic field, corresponding to the allowed values of the total nuclear spin I = 0, 2, 4, 6 and its 
21 + 1 projections Mj for each value of /. The total energy splitting is ~ h x 270 kHz at zero 
field 21 (see the Methods section). Importantly, there is only a single Mj = M F = 6 sublevel 
of \v = 0, J = >, namely the \I = 6, Mj = 6 > level. This level we selectively populate by 
exploiting the dipole selection rule Am f = for linear polarization along the axis of quantization. 
It is the lowest-energy hyperfine sublevel and hence the absolute energy ground state of the Cs 
dimer for magnetic fields above ~ 13 ml 

There are two possibilities for optical transfer from |1 > to |5>. Sequential STIRAP (s- 
STIRAP) uses two consecutive two-photon STIRAP processes, first from |1> to |3> and then 
from |3> to |5>. The second scheme generalizes STIRAP 19 22 to the five-level system 23 : Four- 
photon STIRAP (4p-STIRAP) relies on the existence of a dark state of the form \D> = (fi 2 ^4|l> 
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— fii^ 4 |3 > + ^1^3)5 >)/A with time-dependent Rabi frequencies Qi = fli(t) for lasers L { , 
i — 1,2, 3, 4, and the appropriate normalization function A = A(t). Similar to standard two-photon 
STIRAP, a counter-intuitive pulse sequence rotates the initial state |1> adiabatically into the final 
state, here |5>. For this, L 2 and L 3 couple the three intermediate levels while L 4 and L\ deliver 
time-dependent overlapping pulses with L 4 preceding L\. Fig. [3}) and e show the timings for 
both schemes including the reverse sequence used for detecting the ground-state molecules after a 
certain hold time th- 

We investigate 4p-STIRAP to \v = 0, J = > by interrupting the transfer sequence after a 
given 4p-STIRAP time r and measuring the number of Feshbach molecules, as shown in Fig. [3^. 
The molecules are transferred to 1 5 > in a single step. No molecules in 1 1 > are detected during 
Th as the remaining Feshbach molecules are cleared by L\ at the end of the transfer. When the 
pulse sequence is reversed, a large fraction of the molecules returns to |1 >. Typically, 30% of the 
molecules are recovered after the full double 4p-STIRAP sequence. Assuming equal efficiencies 
for both transfers, the single-pass efficiency is 55%. The inset shows the double 4p-STIRAP 
efficiency versus detuning A 4 of L 4 from the (|4 >— » |5 >)-transition with all other lasers on 
resonance. With ground-state transfer efficiencies between 55% and 60%, about half of the lattice 
sites are occupied by a ground-state molecule. The solid lines in Fig. [3^ represent a simulation of 
4p-STIRAP that takes into account excited-state spontaneous decay and laser linewidth. Transfer 
times are typically 4 fis to 10 /zs. The simulation yields that the transfer efficiency is currently 
limited by a combination of laser linewidth, which is about 10 kHz when averaged over 1 s, and 
imperfect adiabaticity due to finite available laser power to drive the extremely weak transitions of 
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the 5-level scheme . Molecules not transferred to 1 5 > as a result of insufficient phase coherence 
or limited adiabaticity are excited to either 1 2 > or 1 4 > by one of the lasers and are hence pumped 
into a multitude of rovibrational levels, which do not couple to the rovibrational ground state. For 
comparison, the double s-STIRAP efficiency from |3 > to \v — 0, J — > and \v = 0, J = 2 > is 
shown in Fig. [3£ and d, respectively. The solid lines represent a calculation matched to the data for 
standard 3-level STIRAP. With 55%-60%, the total (|1 >-> |5 >)-transfer efficiency for s-STIRAP 
is comparable to 4p-STIRAP. 

A crucial prerequisite for efficient ground-state transfer without heating is perfect matching 
of the motional wave functions for the initial weakly-bound state and the final ground state. A 
mismatch leads to unwanted excitation of higher lattice vibrational levels or bands and hence to 
loss of state control. The lattice thus has to be operated at the magic wavelength condition 25 , i.e. at 
a wavelength that gives equal light shifts for the initial and the final molecular states. We measure 
the lattice band structure and determine the molecular polarizability of the ground- state molecules 
(see the Methods section). Molecules residing in the lowest band of the lattice are excited to the 
first (second) band by phase (amplitude) modulation of the light generating the lattice. Fig. 0] 
shows the measured band energies together with a calculation of the band structure as a function of 
lattice depth. On resonance, excitation to higher bands can readily be observed in momentum space 
as shown in inset a. For comparison, off-resonant modulation transfers hardly any population into 
higher bands (see inset b). We determine the band energies by taking modulation spectra as shown 
in inset c. We then use the band structure calculation to fit all measured resonance positions with 
the molecular dynamical polarizability P\ v =o> as the single free parameter. These measurements 
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aredonefor |f = 0, J = 2>. We obtain P\ v =o > = 2.1(1) x P a , where P a is the atomic polarizability. 
For the initial, weakly-bound Feshbach molecules in level \g > we obtain P\ g> = 2.0(1) x P a . 
Hence the magic wavelength condition is well fulfilled. 

We measure the lifetime r of the molecules in the lattice by varying the hold time r h for up 
to 20 s and recording the number of remaining molecules as shown in Fig. [5j To reduce inelastic 
light scattering, the lattice depth was adiabatically reduced to about 41.5 Er after the 4p-STIRAP 
transfer. An exponential fit gives a lifetime of r = 8.1(6) s. We attribute this long lifetime 
to the large detuning A L 6.9 THz from the lowest 0+ level with predominant A X Y^ singlet 
contribution as shown in the inset to Fig. [5J Levels of the 0+ system that lie below this are almost 
purely of b 3 U u character and thus make negligible contributions to the optical excitation rate. 

We are now in a position to determine collisional properties of ultracold ground-state molecules 
in a fully state-selective way. At magnetic fields beyond 13 mT, where the level \I = 6, Mj = 6 > 
becomes the absolute ground state, the sample should show collisional stability and thus allow 
the formation of a BEC of ground-state molecules when the lattice is adiabatically removed 12 . 
The long coherence times and the perfect decoupling from the environment in an optical lattice as 
demonstrated here will enable a new generation of precision measurements 26 ' 27 . Furthermore, our 
results can readily be generalized to heteronuclear systems such as KRb 14 and RbCs 28 , opening up 
the possibility to study dipolar quantum phases in optical lattices. 
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1 Methods 



Lattice loading We first follow the procedure detailed in Ref. 18 . In brief, we produce an atomic 
BEC with typically 1 x 10 5 Cs atoms in the lowest hyperfine sublevel F a = 3, rriFa = 3 in a 
crossed optical dipole trap. We then drive the superfluid-to-Mott-insulator phase transition 5 by 
exponentially ramping up the power in a three-dimensional optical lattice within about 300 ms. 
The lattice is generated by three mutually-orthogonal, retro-reflected laser beams at a wavelength 
of A = 1064.5 nm, each with a 1/e-waist of about 350 /mi. While ramping up the lattice potential, 
the power in the two dipole-trap beams is increased to ensure that the central density in the trap is 
sufficiently high to allow formation of atom pairs at the central wells of the lattice, but not too high 
to lead to triply occupied sites. We ramp the lattice to a depth of about 20 E R before Feshbach 
association. Here, E R = h 2 / (2m a \ 2 ) = fc B x 64 nK is the atomic photon-recoil energy with the 
mass m a of the Cs atom and Boltzmann's constant k B . Up to 45% of the atoms reside at doubly 
occupied lattice sites. We estimate this number from the molecule production efficiency. This 
value is close to the maximum of 53% that can be achieved for loading a harmonically confined 
BEC in the Thomas-Fermi limit into an optical lattice 6 . 

For the molecules, the recoil energy is E R = h 2 / (2 x 2m a \ 2 ). The polarizability of Feshbach 
molecules is twice the atomic polarizability. The same lattice that has a depth of 20 E R for the 
atoms has thus a depth of 80 E R for the Feshbach molecules. 



Feshbach association and Feshbach state transfer We efficiently produce weakly bound Cs 2 
Feshbach molecules in the presence of the optical lattice by a magnetic field sweep 7 across a 
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narrow g-wave Feshbach resonance with its pole at a magnetic field value of B = 1.98 mT '■ . 
The molecules are initially in level \g >, for which £ = 4. Here, £ is the quantum number associated 
with the mechanical rotation of the nuclei 7 . We subsequently transfer the molecules via level \g 2 > 
with 95% efficiency into level \ s>= 1 1 > with £ = by magnetic field ramping 17 ' 18 as shown in Fig. 
[2J5. For this level, the transitions to excited molecular levels are stronger than for the initial level 
\g > 30 . We obtain up to 2.5 xlO 4 Feshbach molecules in the lattice in the desired starting state. 
We detect the molecules in 1 1 > by reversing the Feshbach state transfer sequence, dissociating 
the molecules at the Feshbach resonance and detecting the resulting atoms by standard absorption 
imaging 29 . 

Molecular states for ground state transfer The relevant molecular states for Cs 2 are shown in 
Fig. [2^. Levels \2 > and |4 > belong to the coupled (A 1 E+ - 6 3 I1 U )0+ potentials 13 . We have 
recently identified suitable transitions linking |1 > to |5>, where levels |2 >, |3 >, and |4 > were 
chosen to give balanced transition strengths on the four optical transitions 24 - 30 . For |3 > we choose 
either \v = 73, J = 2 > or \v = 73, J = > of the ground state with a binding energy of 

~ hex 1061 cm -1 . 

Hyperfine structure of the rovibronic ground state The hyperfine levels are calculated using 
the molecular constants from Ref. 21 by constructing and diagonalizing a Hamiltonian matrix in 
an uncoupled basis set of functions representing the molecular rotation and the spins of the two 
nuclei, using the matrix elements given in the Appendix of Ref. 21 . For J = states the hyperfine 
structure is dominated by the scalar spin- spin coupling and the nuclear Zeeman effect, but for 
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J > additional terms are important. 

STIRAP laser setup STIRAP lasers L { with % — 1, 2, 3, 4 are continuous-wave grating- stabilized 
tunable diode lasers, which are stabilized to optical resonators for short-term stability and refer- 
enced to an optical frequency comb for long-term stability and reproducibility. We estimate the 
linewidth of the lasers to be about 10 kHz. In order to ensure minimum momentum recoil imparted 
on the molecules, the beams for lasers L\ and L 2 are co-propagating. The beams for L 3 and L 4 are 
also co-propagating but run antiparallel to the beams of L 1 and L 2 . All beams run horizontally and 
are linearly polarized with the polarization axis in the vertical direction, parallel to the direction of 
the magnetic field, which defines the axis of quantization. We operate at Rabi frequencies in the 
range of 2irx (1 to 4) MHz. 

Polarizability measurement For determining the ground state molecular polarizability, transfer 
to \v = > is performed at a fixed lattice depth of 83 E R for \v = > molecules. The lattice 
depth is then ramped to the desired value within 50 ms. For phase modulation of the lattice, the 
frequency of one lattice beam is usually modulated with a modulation depth of 2 MHz at the 
desired frequency for about 10 ms. For amplitude modulation, the intensity is typically modulated 
by 20% for about 10 ms. 
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Figure 1 : Molecular quantum gas preparation procedure. A BEC of Cs atoms is loaded 
into an optical lattice. By increasing the lattice depth, a Mott-insulator state (MI) with preferen- 
tially two atoms per site is created. Feshbach association (FA) subsequently converts atom pairs 
into weakly-bound molecules in state |1>. These are then transferred in the presence of the lat- 
tice to a specific hyperfine level \I = 6, Mj = 6> of the rovibronic ground state |5> = X l Y>+ 
\v = 0, J = > by a stimulated four-photon process (STIRAP) involving lasers L { with Rabi 
frequencies fli % — 1, 2, 3, 4, and three intermediate levels |2>, |3>, and |4>. 
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Figure 2: Molecular potentials and level schemes for ground-state transfer, a, The four- 
photon transfer from the weakly-bound Feshbach level |1>= \v 155 > (not resolved near the 
6S+6S asymptote) to the rovibrational ground state |5>= \v = 0, J = 0> of the singlet X 1 ^ 
potential involves the deeply bound level |3>= \v = 73 > of the X 1 ^ potential 13 and the levels 
|2>= |i/ = 225, J = l> and |4>= |i/ = 61, J = l> of the electronically excited - 6 3 n„) 0+ 

potentials 24 30 . The laser wavelengths for L x , L 2 , L 3 , and L 4 are near 1126 nm, 1006 nm, 1351 
nm, and 1003 nm, respectively, b, Zeeman diagram for weakly bound molecules near the 6S+6S 
asymptote. Molecules are associated at a g-wave Feshbach resonance 29 at 1 .98 mT (FA) and then 
transferred to the desired starting level |1 >= \s > for optical transfer via three avoided level 
crossings by slow (arrows 1,3) and fast (arrow 2) magnetic field ramps 17 . The binding energy is 
given with respect to the (F ai = 3, m Fai = 3) x (F a2 = 3, m Fa2 = 3) two-atom lowest hyperfine 
asymptote. All Feshbach levels are characterized by M F = 6. c, Calculated Zeeman diagram for 
the hyperfine manifold of the rovibronic ground state |5>= \v — 0, J = 0>. The optical transfer 
goes selectively to level |7 = 6, Mj = 6>, indicated in red. This level becomes the lowest-energy 
absolute ground state for magnetic-field values above ~ 13 mT. There are no avoided crossings 
between different hyperfine sublevels 21 . 
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Figure 3: STIRAP transfer to the rovibronic ground state |5>= \u = 0, J = 0> and 
back, a, 4p-STIRAP transfer and b, schematic timing for the Rabi frequencies Q i9 % — 1,2, 3, 4: 
Number of molecules in state 1 1> as a function of 4p-STIRAP time r for all 4 lasers on resonance. 
The lattice depth is 80 E R and 83 E R for molecules in levels |1> and |5>, respectively. Data 
points represent a single experimental realization, not an average over several runs. The solid 
line is a 4p-STIRAP model calculation. r h is the hold time in |5 >= \v — 0, J = >. Inset: 
double 4p-STIRAP efficiency as a function of the detuning A 4 of laser L 4 and corresponding 
model calculation. The peak corresponds to a single-pass efficiency of 57%. c, and d, s-STIRAP: 
Double STIRAP efficiency for the inner two-photon STIRAP from |3 > to \v = 0, J = > (c) and 
to \v — 0, J = 2 > (d) and back, corresponding to the dotted bar in the timing sequence in e, as 
a function of the detuning A 4 of laser L 4 . The number of molecules is normalized to the initial 
number in |3 >. All measurements are performed at an offset magnetic field of 1.9 ml 
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Figure 4: Lattice band structure for \v — > molecules. Band energies as a function of 
lattice depth in units of the molecular recoil energy E R as measured by phase and amplitude mod- 
ulation of the lattice. The lattice bands are labeled by (k, I, m), where k, I, and m give the number 
of vibrational quanta along the three spatial directions in the limit of a deep lattice. The horizon- 
tal position of the data points (filled circles, representing the position of excitation resonances as 
shown in inset c) is given by the molecular polarizability, which is determined by a fit of the data 
to the band structure. Inset a shows the molecular momentum distribution after transfer to higher 
lattice bands by resonant lattice amplitude modulation. The distribution represents an average of 5 
experimental runs, smoothed with a Gaussian filter. For comparison, inset b shows that hardly any 
population is transferred to higher bands for off-resonant modulation. The rectangular shape of 
the first Brillouin zone is easily seen. Inset c shows typical excitation spectra for amplitude (top) 
and phase (bottom) modulation at 83 E R . For these, we determine the number of molecules in the 
first Brillouin zone as a function of the excitation frequency. The solid lines are Gaussian fits. The 
resonance at 22.5 kHz corresponds to excitation to the nearly-degenerate bands (2,0,0) and (1,1,0) 
(not resolved). The resonance at 20.6 kHz is a two-phonon excitation to (4,0,0). 
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Figure 5: Lifetime of trapped ground-state molecules in the optical lattice. Normalized 
number of molecules in|5>= |z/ = 0,J = 0>asa function of hold time r^. The solid line is 
an exponential fit, yielding a lifetime of 8.1(6)s. The inset schematically shows the excited-state 
potentials to which off-resonant optical excitation is possible (cf. Fig. [2^l). A l is the detuning of 
the lattice light at 1064.5 nm with respect to the lowest 0+ level with A 1 ^ character. During the 
hold time, all STIRAP laser fields are turned off. 
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